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Abstract 

A family of sequences produced by a non-homogeneous linear recurrence formula 
derived from the geometry of circles inscribed in lenses is introduced and studied. 
Mysterious "underground" sequences underlying them are discovered in this paper. 
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1 Introduction 

We investigate a new family of integer sequences. They are generated by a geometric con- 
struction, which we now describe. 

Start with three circles of curvatures a, b, and c centered on the same line, so that pairs of 
consecutive circles are tangent, as in the left side of Figure 1.1. The three circles determine 
a pair of congruent circles that are simultaneously tangent to the original triple (right side 
of Figure 1.1). The common region formed by this pair defines a symmetric lens. Now, 
continue to inscribe circles inside the lens, as shown in Fig. 1.1. The resulting chain of 
circles defines a bilateral sequence of curvatures (foj), i E Z. Sequences obtained this way 
will be called lens sequences. 



^ bo=a bi=b b2=c bs ^4 bs 



Figure 1.1: Three circles determine a sequence 

Terminology. The resulting circles form a lens circle chain. We shall say that a triplet 
of circles (a, b, c) generates the lens sequence, and we will call it a seed of the sequence. 
Notice that any three consecutive terms of a lens sequence define a seed. The two circles 
that form the lens will be called lens circles. 



Notation. Typically, we denote circles and their curvatures by the same letter (circle a 
has curvature a, i.e., radius 1/a). 

Our opening result (proved in the next section) is this: 

Theorem A. Let a, b and c be the curvatures of the initial three circles generating a lens 
sequence, b ^0. Then the sequence is determined by the following inhomogeneous three-term 
recurrence formula: 

bn = (ybn-l-bn-2+ P , (1-1) 

where a and (3 are constants determined by the original triple: 

ab + bc + ca , ^ b"^ — ac , , 

a = 1 and p= — . 1.2 

0^ b 

In particular, if bo = a and bi = b then 62 = c . 

Constants a and P are "invariants" of the sequence - their values may be determined 
from any three consecutive terms of the sequence. 

If a, b, c as well as a and /? are integers, then (6„) is an integer sequence. Surprisingly, 
this family of sequences includes a wide range of known sequences [9]. However, for some of 
these sequences, properties that we develop in this paper seem to be new. For now, let us 
look at a few examples. 

Example 1.1 (Vesica Piscis). Starting with (a, b, c) = (3,1,3) we get the recurrence formula 

bn = 146„„i - bn-2 - 8 , 

which produces 

... 3, 1, 3, 33, 451, 6273, 87363, 1216801, 16947843, 236052993, . . . , 
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a sequence known as A011922 . Note that if one starts with values twice as large, the whole 
sequence is doubled: 

. . .6, 2, 6, 66, 902, 12546, 174726, 2433602, 33895686, 472105986, . . . 

and the recurrence formula becomes 6„ = 146„_i — bn-2 — 16 (the same a and twice f3). (The 
lens circles traverse each others center, forming a well-known figure of Vesica Piscis, hence 
the name of the example.) 

Example 1.2 (Golden Vesica). Start with (a, 6, c) = (1,2,10). Equations (4.2) give a = 7 
and /? = — 3; hence the sequence is generated by 

bfi = 76n-i ~ bn-2 ~ 3 

and is 

1, 2, 10, 65, 442, 3026, 20737, 142130, 974170, . . . 

for positive n. This sequence, listed as Sloan's A064170 , is known for its interesting prop- 
erties. Its terms are products of pairs of non-consecutive Fibonacci numbers: 1-2, 2-5, 
5 • 13, 13 • 34, . . . , etc. They also coincide with the denominators in a system of Egyp- 
tian fraction for ratios of consecutive Fibonacci numbers: 1/2 = 1/2, 3/5 = 1/2 -|- 1/10, 
8/13 = 1/2 -|- 1/10 + 1/65, etc. (The geometry of the lens relates to the golden proportion, 
hence the proposed name.) 

Example 1.3. Triplet (—1,3,15) gives 

... 99, 63, 35, 15, 3, -1, 3, 15, 35, 63, 99, 143, . . . 

from the recurrence 

bn = 26n-l — bn~2 + 8 . 

The sequence (3, 15, 35, . . . ) is known as A000466 and is defined by 6„ = 4n^ — 1. The 
occurrence of negative curvatures will be explained later. 

Example 1.4. A lens sequence does not necessarily need to be symmetric. For instance the 
triple (2,1,3) produces the following bilateral sequence: 

. . . , 12972, 1311, 133, 14, 2, 1, 3, 24, 232, 2291, 22673, 224434, . . . 

Example 1.5 (More Examples). The lens sequences possess an ample diversity. They 
include such basic examples as (i) the powers of 2 ( A000079 ), and (ii) triangular numbers 
(A000217). 

(i) 1, 2, 4, 8, 16, 32, 64, 128, 256, ... bn = 5/26„_i - 6„_2 

(ii) 1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 66, 78, 91, . . . 6„ = 26„_i - bn-2 + 1 

A more extensive list with references to OEIS [9] is provided in Tables 4.1 through 4.3. 
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Figure 1.2: Sequences A011922 , A064170 , and A000466 in the Apollonian window 



Remark 1.6. Certain circle packings, known as Apollonian gaskets [5], result in integral 
curvatures for all of the circles (see e.g., [4]). One such gasket, an Apollonian Window (see 
[2, 3]), is presented in Fig. 1.2. Interestingly, it contains an infinite number of lens sequences, 
from which the three shown in Fig. 1.2 are especially conspicuous. They correspond to the 
Examples 1.1, 1.2, and 1.3, given above (up to scaling). This observation was the author's 
initial motivation for this study. 

The negative term of Example 1.3 (the third one in Figure 1.2) has a clear geometric 
meaning: it is the curvature of the greatest disc in the sequence, which — unlike the other 
discs in the sequence — contains the lens circles (as well as the rest of the sequence). 

A number of interesting features are common to all lens sequences: 



1. Limits. In many cases, the limit of the ratios of consecutive entries is well defined. For 
instance, referring to the above examples: 

Example 1 [ API 1922] : lim ^ = 7 + 4^3 = (2 + ^3)^ 

n^oo On 



Example 2 [ A064170] : lim 



n^co bn 2 



b 

Example 3 [ A000466 ]: lim = 1 

n^oo bn 



These numbers are examples of Pisot numbers and will be called characteristic constants of 
the sequences denoted A = (a + Wa"^ — 4). 



2. Sums. The reciprocals of curvatures are the circles' radii. Their sum is determined by 
the length of the lens. For instance: 
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Example 1 [ API 1922 ]: E^i l/k = 1 + 1/3 + 1/33 + 1/451 + . . . = ^3/2 

Example 2 [ A064170] : l/k = 1 + 1/2 + 1/10 + 1/65 + . . . = (1 + ^5)/2 = 1.618 . 

Example 3 [ A000466] : S^^^ l/6„ = J^n - 1) = V3 + 1/15 + 1/35 + . . . = 1 

3. Binet-type formulas. For a 7^ 2, the curvatures may be expressed in terms a non- 
homogeneous Binet-type formulas: 

Example 1 [Vesica Piscis, A011922 ]: 3 1 3 33 451 6273 87363 . . . , 

_ 4 + (2 + v^)^" + (2 - VS)^^ 

bn- g 

Example 2. [Golden Vesica [ A064170] : ... 2 1 2 10 35 442 ... , 



3+(^)°+( 



in / _\ 4n 



Example 3. Also non-symmetric lens sequences can be expressed this way. For instance, the 
sequence extended from (6, 2, 3), which is (. . . 2346 299 39 6 2 3 15 110 858 6747. . . ), with 
the recurrence 6„ = 86„__i — 6„_2 — 7, may be obtained from 



^ _ (25 - 3V15)(4 + yiS)" + (25 + 3V15)(4 - ^15)" 7 
" ~ 60 ^ 6 

The above properties are known for some of the sequences, but now they acquire a 
geometric interpretation. Other related concepts include geometry of inversions, Chebyshev 
polynomials, etc. 

The most remarkable and perhaps surprising property is that the integer lens sequences 
are "shifted squares" of yet deeper integer "underground" sequences. This discovery is the 
topic of the final section of this paper. 



2 Recurrence formula from geometry 

In this section we prove Theorem A on the recurrence formula for lens sequences. A reader 
interested in the algebraic properties of these sequences may skip it without loss of continuity. 

We shall need a theorem on circle configurations generalizing that of Descartes' theorem 
on "kissing circles" ([1, 2]). If Ci and C2 denote two circles of radii ri and r2 respectively, 
and d denotes the distance between their centers, then one defines a product of the circles 
as 
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which we propose to call the Pedoe product. Its values for a few cases are shown in Figure 
2.1. For any four circles Ci, i = 1, ... ,4, define a configuration matrix / as the matrix 
with entries 

fij = {Ci, Cj) , 

where the brackets denote the Pedoe inner product of circles. 



intersecting 
circles 




external 
tangency 



9 = 0° 




{Cj,C2) = COS 9 



{Cj,C2) = 1 



internal 
tangency 




9= 180° 



{Cj,C2) = -1 



Figure 2.1: Pedoe product of circles 

Theorem 2.1 ([2]). A configuration of four circles in general position satisfies the following 
quadratic equation 

b^Fb = (2.2) 

where b = 62, &3, ^4]"^ is the vector made of the curvatures of the four circles, and where 
F = f^^ is the inverse of the configuration matrix. 

Remark 2.2. Equation (2.2) is only a fragment of the full matrix formula, which incorpo- 
rates also the positions of the centers of the circles. For more on this theorem, its proof, 
and the associated Minkowski geometry of circles, see [2]. For our purposes Formula (2.2) is 
sufficient. 

Notation. In the following by a lens we mean "symmetric lens" — the intersection of the 
interiors (exteriors) of two congruent circles, called in this context lens circles. A chain of 
circles is a sequence of circles such that every two consecutive circles are tangent. 

We are now ready to prove the basic result. 

Theorem 2.3. A sequence (bn) of curvatures of a chain of circles inscribed in a lens satisfies 
a non-homogeneous linear recurrence formula of the form 

bn+i = abn -bn-i+ P 

for some constants a and [3, with 

where K is the Pedoe product of the two lens circles and A = 1/R is the curvature of each 
lens circle. 
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Proof. Consider two consecutive circles in the lens, of curvatures say a and b. Denote the 
curvatures of the circles that form the lens by A, and their Pedoe product by K [K = cos ip, 
if the circles intersect). The configuration matrix / and its inverse are easy to find. In the 
case of converging lenses we can read it off from Fig. 2.2a: 





■ -1 K 


-1 


-1 " 


/ = 


K -1 


-1 


-1 


-1 -1 


-1 


+ 1 




-1 -1 


+ 1 


-1 



where the indices are ordered as (A, A, x, y). Its inverse F is easy to find, and the master 
equation (2.1), after multiplying by a factor of 8, becomes: 



A 
A 

X 

y 



n T 



K+1 K+1 

-4 4 



2 
2 



K+1 K+1 

2 2 K + 1 K 



A 
A 

X 

y 



2 2 K-?> K + l 
This quadratic equation is equivalent to: 

(1 + K)x'^ + (1 + K)y'^ + 2{K - ?,)xy + SAx + SAy = . 



One may solve it for y to get two solutions (corresponding to two signs at the square root): 

AA + iK- 3)a; ± 2^/2(1 - K^x"^ - SAx + 4^2 ^ ^ 

1/1.2 = ^ (2.4) 

Note that the two solutions yi and ?/2 correspond to the two possible circles tangent to x: 
one on the left and one on the right. To eliminate radicals, add the two solutions: 



yi + y2 



6-2K 



8A 



-X — 



l + K l + K 

Since the triple {yi,x, 1/2) forms a sequence in a chain of inscribed circles, we may label these 
curvatures as = yi, hn = x, and bn+i = 2/2, to get 

bn+i + hn-i = abn + P, 

which is equivalent to (2.3). The case of the diverging lens results from similar reasoning, 
with slightly different initial matrix F. □ 

Corollary 2.4. The sequence constants are related: a + Rj3 = —2. 

Now let us see how three circles determine a sequence. 

Theorem 2.5. Let a,b and c be curvatures of three consecutive circles inscribed in a lens. 
Then the sequence of the circle curvatures is determined by the following three-term recur- 
rence formula: 

bn = abn-i-bn-2 + P , (2.5) 



where a 



ab+bc+ca 



■1 and (3 



f,2- 



//6o = 0- andbi = b then 62 = c, that is, ab+f3 = a+c. 



b2 b 

Proof. We apply Theorem 2.1 in each of the three steps to a different quadruple of circles. 
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Figure 2.2: The three steps of the proof of Thm. 2.5 



Step 1. Consider a configuration of four circles: a triple of three consecutive circles in the 
chain, say a, b, c, plus one circle forming the lens, say d = 1/R. The Pedoe product of two 
external circles a and c may be easily evaluated; since the distance between their centers is 
(1/a + 2/6 + 1/c), we have: 



(a,c) 



ii+i+',f-iif-iiy 



2 ■ 



1 1 



ab + be + ca 
b^ 



Denote the main fraction of the last expression by 2; = "''+^^+^" . Then the configuration 
matrix and its inverse are 





-1 


1 


2z + l 


-1 " 


/ = 


1 

2z + l 


-1 
1 


1 

-1 


-1 
-1 




-1 


-1 


-1 


-1 



1 



]_ 

z+1 

1 -{z + 1) 
-1 z-1 



1 

z+1 
1 

]_ 

z+1 

-1 -(^ + 1) 



-1 

z-1 

-1 



(the order of entries is: abed). Denoting v = [a,b,e,d]'^ and solving the quadratic equation 
v'^Fv = for (i readily leads to 



d 



b{ae — fe^ 



ab + be + ea + 6^ 



This gives us the curvature of each of the two lens circles. Now we need to find the product 
of these two lens circles. 

Step 2. Use a quadruple of circles: a, b, and the two circles forming the lens, d and d'. The 
latter two have the same curvature d = d', the value of which we know from the previous 
step. The goal is to find the Pedoe product K = {d, d') . The configuration matrix and its 
inverse are 





-1 


1 


-1 


-1 


/ = 


1 


-1 


-1 


-1 


-1 


-1 


-1 


K 




-1 


-1 


K 


-1 



-1-K 
3-K 

-2 

-2 



i-K 
-1-K 

-2 

-2 



-2 
-2 



K+l 
4 



K+l 



K+l 
4 

K+l 
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(the order of indices agrees with abdd'). Applying vector v = [a,h,d,d\^ to the quadratic 
equation v"^Fv = gives 

01,2 

K = — ^ - 1 (2.6) 

Step 3. Now we can either build the matrix for configuration (c) in Figure 2.2 and mimic 
the proof of Theorem 2.3, or simply substitute for K from (2.6) in (2.4) to get the result. □ 

3 More on lens geometry 

Although we are mainly interested in the algebraic properties of lens sequences, some geo- 
metric properties explicate their algebraic behavior. Below, we summarize basic facts. 

Proposition 3.1. The radius R of the lens circles is determined by three circles and may 
be expressed in terms of the sequence constants a and /3: 

a + 2 (a + &)(& + c) 

The Pedoe inner product of the lens circles is 

6 - a 1 / 5 \ ^ , , , 

K=- = - — --1 = - - -l = cos¥P, (3.1b) 



2 + a (a + 6)(6 + c) 2 

where the last equation is valid if the circles intersect. The length L of the lens, if defined, 
is 



^_op _ 2v/^^ _^ v/(a + &)(& + c)[(a + &)(& + c)-4&2] 

Va + 2 /3 (ac-62)6 ' ^ ' ^' 

The separation of the lens circles (distance between their centers) is 

v/aT2 /3 R y^T2 

Proof. All are direct corollaries of Theorem 2.3 and simple geometric constructions. □ 
Figure 3.1 contains these findings for easy reference. 

Figure 3.2 categorizes a variety of geometric situations for a lens sequence. In the case of 
converging lenses, when two circles of radius R intersect at angle ip, the recurrence formula 
is: 

- 2 I bn-l - bn-2 - 



1 + COS ip J Rl + COS (p 

This answers the question of which lenses may lead to integer sequences. Indeed, denote 
n = z—^ — • Then a = n — 2, 3 = —n/R. For n to be an integer, n G N, we need 
COS if = 8/n — 1. Table 3.1 shows some values. 
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Figure 3.1: Sequence constants and geometry of a lens 

converging lens diverging lens diverging lens prism slab 

(regular) (corrupted) 

a>2, p<0 a<2, p>0 a>2, p>0 a>2, p>0 a=2, p=0 

Figure 3.2: Types of lenses and associated sequences 

Figure 3.2 relates the geometry of lenses to the values of the sequence constant a and 
the Pedoe product K. Inspect Tables 4.1 through 4.3 (and associated figures) in the next 
section for various examples of lens sequences. 

Note that if a ^ 2, then only external sequences are possible (corresponding to diverging 
lenses). Moreover, if a < 2, then the integer sequence must be periodic. If a > 2, then we 
can have two families of sequences: inner (inside a converging lens) or outer (outside the 
lens circles, i.e., inside a corrupted diverging lens [corrupted, because of the missing central 
part)). In the case of the outer sequence we will have exactly one negative entry (the most 
external circle) or two adjacent "0" entries (two vertical lines). 
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n 


K 


a 




example 


n 
u 


OO 


9 


pcllOQlC I^OICICI ZJ 




1 


7 


-1 


periodic (order 3) 


(6,3,-2) 


2 


3 





periodic (order 4) 


(3,6,2,-1) 


3 


5/3 


1 


periodic (order 6) 


(2,10,15,12,4,-1) 


4 


1 


2 


Example 1.3 


(-1,3,15) 


5 


3/4 


3 




(3,2,2,3) 


6 


1/3 


4 




(2,1,1,2) 


7 


1/7 


5 




(5,2,2,5) 


8 





6 


orthogonal 


(3,1,1,3) 


9 


-1/9 


7 


golden Vesica 


(2,1,2) 


10 


-1/5 


8 




(4,1,1,4) 


11 


-3/11 


9 






12 


-1/3 


10 




(3,1,2) 


13 


-5/3 


11 






14 


-3/7 


12 




(6,1,1,6) 


15 


-7/15 


13 




(4,1,2) 


ID 


1 /9 


1 zL 


V cbiCd Jr ibCib 




17 


-9/17 


15 






18 


-5/9 


16 




(8,1,1,8) 


20 


-3/5 


17 







disjoint lens 
circles 

tangent lens 
circles 



intersecting lens 
circles 



Table 3.1: Admissible values of a 

4 Basic algebraic properties of lens sequences 

Theorem A suggests the following definition: 

Definition 4.1. A formal sequence extended from a triplet (a, 6, c), called a seed, is 
defined by the following inhomogeneous three-term recurrence formula: 



hn^a hn-l - bn-2 + P , 

where a and /? are constants determined by the original triple: 

ab + bc + ca 



a 



— 1 and P 



6^ — ac 



(4.1) 



(4.2) 



and bo — a and bi — b. (It follows that b2 — c) . 

The values of the constants a and P do not depend on the particular choice of the triplet 
of consecutive terms (seed). Moreover, if a > —2, then the sequence may be interpreted in 
terms of a chain of circles inscribed in a lens made by two disks each of curvature — 
—(3 /[a + 2) separated by distance 5 — AR/^/a + 2. Formal integer lens sequences exist also 
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for a < —2 (see Table 4.4 for examples), but in such a case the geometric interpretation is 
unclear as the distance between the lens circles becomes imaginary. 

In general, lens sequences take real values. However, if any three terms of a lens sequence 
are rational, so is the whole sequence. The question whether a particular seed produces an 
integer sequence will be answered for now this way: 

Proposition 4.2 (Integrality Criterion 1). If b\ac and b'^\{ab + bc + ca), for any a,b,c & N, 
then a lens sequence extended from (a, b, c) consists of integers. 

Here are basic properties of lens sequences: 

Proposition 4.3. Let (6„) be a lens sequence. Then the following holds: 

(i) Sequence {bn) satisfies a homogeneous 4-term linear recurrence formula 

bn = {a + l)6„_i - (a + 1)6„_2 + bn-3 ■ (4.3) 

(ii) If a > 2 then the sum of the reciprocals converges and equals: 



71= — OO 



(iii) If a >2 then the limit of the ratios of consecutive terms exists and equals: 

a + - 4 ^/a + 2 + ^/a-2 
A = ^ = ^ . (4.5) 

(iv) If a 7^ 2, then the lens sequence generated from a seed (a, 6, c) has the following 
Binet-like formula 

bn = wA" + wA" + 7 (4.6) 



where 



and where 



a + y/a'^ — 4: - a ~ _ 4 

A = A = . 

2 2 

a — 2b + c c — a 1-^ — /3 

w = — ; ^ — I — rV a — 4, 7 = 

2(a-2) 2{a^-4) ' ' a -2 



and w and A denote conjugates of w and A in Q(va^ — 4), respectively. In particular, 
(a,6,c) = (6_i,6o,&i)- 

Proof, (i) Elementary, (ii) From the geometry of lenses, cf. (3.1c). See also Figure 4.1. 
(iii) Divide the recurrence formula by 6„_i to get 

bn/bn-i = a~ bn-2/bn-l + P/bn-1 ■ 

For large values of n, since the sequence is divergent, the last term becomes irrelevant and 
the equation becomes A = a — 1/ A, or simply 

A2-aA + l = 0, (4.7) 
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Figure 4.1: Geometric meaning of ^ l/6„ and A (see Proposition 4.3) 

n 



with the solution as above. Figure 4.1 provides the geometric insight, which also relates A to 
the lens angle via similar triangles, (iv) Define a new sequence whose entries are shifted by a 
constant, namely a„ = 6„ + /?/(« — 2). The sequence (a„) satisfies a homogeneous three-term 
recurrence formula which resolves to 4.6 by the standard procedure. □ 

The value of A (given by (4.5)) will be called the characteristic constant of the se- 
quence. The ring over rational numbers generated by a/^^ — 4 plays an important role in 
other properties of lens sequences, as we shall soon see. Note that the sequence constant a 
may be expressed in terms of the characteristic constant in a graceful way: 

a = A + — 



Alternative generating formulae 

Note that the three term formula (4.1), with given coefficients a and (3, requires only two 
initial entries to produce a sequence. Yet not all such initial values will produce a lens 
sequence of the type under discussion. This is because arbitrary initial values bo and bi do 
not need to be geometrically inscribable into a lens defined by a and (3 as two consecutive 
circles. The following will clarify the situation: 

Proposition 4.4 (Compatibility condition). Two consecutive circles a and 6 in a lens chain 
satisfy the following condition: 

0^ + 52 = _^ _^ 5) _ (4 8) 

Proof. Eliminate c from the expressions for a and (3 in (4.2), and simplify. □ 

The above formula may actually be used as an alternative definition of lens sequences. 
Indeed, 

Proposition 4.5. Consider the following properties: 

(a) Recurrence bn+i = abn — bn-i + P , for all n, 
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(b) Constants a = f'r.-ibn+b„b^n+i+b„+ibr.-i _ ^ ^ _ bl-b^^_^K+^^ 

(c) a"^ + b"^ = aab + (3{a + 6) , for a = 6„ and h = 6„+i. 
The following descriptions of a sequence (bi) are equivalent: 

(i) Recurrence (a), and constants (b) for some n (definition of a lens sequence); 

(ii) Recurrence (a) and compatibility condition (c) for some n; 

(iii) Any of the two constant formulas (b) for all n; 

(iv) Compatibility condition (c) for all n. 

Proof. Let us show that the compatibility condition together with the recurrence theorem 
imply our standard formulas for both a and /3. Starting with (ii) we get 

bl+i + bl = a bnbn+i + /3(6„ + 6„+i) 
= bn+i{abn + P) + Pbn 

= bn+l{bn+l + bn-l) + . 

Subtracting b'^j^^ from both sides, we get 6^ = bn+ibn-i + Pbn, which gives 

^ _ b\ ~ bn-lbn+l 

K 

To get a, substitute this result in the compatibility condition and simplify. The other 
equivalences follow easily. □ 

Here is yet another intriguing formula that will prove itself handy later. 

Proposition 4.6. Constant a has an alternative form involving any four consecutive entries 
of a lens sequence: 

« = — + — . (4.9) 
Proof. Start with the formula for f3 and express it as follows: 

P = 7 = bn 7 =^ bn- 13= . 

bn bn bn 

Use the recurrence formula to modify the left-hand side of the last equation, 

, , _ bn+lbn-l 

bn 

Now extract a to get 4.9. □ 
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Integrality condition 



Lens sequences are self-generating in the sense that any three consecutive entries (a,6, c), 
a seed, determine the whole sequence (unless 6 = 0). The question is how to choose seeds 
(a, b, c) in order to obtain lens sequences that are integer. Below we give only a partial 
answer to this problem of integrality conditions; the last section will provide the solution. 

Definition 4.7. An integer lens sequence is primitive if the common divisor of three con- 
secutive entries is 1. 

Proposition 4.8. If gcd(6fc, &fc+i, bk+2) = n holds for some k, then it holds for all A; G Z. 

Proof. If n divides each of (6^, &fc+i, bk+2), then it divides /? of the recurrence formula. Hence 
it divides the neighboring terms 6^+2 and By induction, n divides every term of the 

sequence. □ 

Recall that to insure that a and (3 are integers, we need to choose (a, 6, c) so that b\ac 
and b'^\ab + be + ca, or, equivalently, that 6^|(a + b){b + c). Thus triples of the form (a, 1, c) 
always generate integer sequences for any a, c G N . For now, let us review the following 
families of integer lens sequences: 

A. A sequence is called central if it contains a triple of the form (a, 6, a). Only b = ±1 
leads to primitive integer sequences. 



B. A lens sequence is called bicentral, if it contains a quadruplet of the form {a,b,b,a). 
Only if b is chosen from {0,1,2}, does a primitive integer sequence result. (For the case 
6 = 0, the seed needs to be chosen in the form (0, a, c)). 

In either case A or case B, the sequence is called symmetric. Table 4.1 shows examples 
of symmetric lenses for small values of the initial terms. (Only the right tail is displayed.) 



C. Here is a method of getting a not-necessarily symmetric integer sequence: choose arbi- 
trarily a couple (a, b) and some integer k (only one of a or 6 can be negative). Then a seed 
(a, b, c) with 

c = b{bk-l), (4.10) 

will generate an integer sequence. Indeed, calculate the recurrence constants a and (3 from 
(4.10) to get the integer values 

a = (a + b)k-2 , , 

(4.11) 

f3 = [a + b) — abk . 

The triple of integers [a, b; k] will be called the label of a lens sequence of this type. (Note 
that it includes the symmetric lens sequences as a special case; for examples of non-symmetric 
sequences, see Table 4.2.) 

The pair (a, b) in the label may be chosen so that it contains the smallest element of the 
sequence. 
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central 

elements constants 



OEIS nr sequence 



label and symbol 



Central sequences 



1. 


(2,1,2) 


a 


= 7,/? = 


-3 


A064170 


1, 2, 10, 65, 442, 3026, 20737, 142130, . . . 


[1,2;3] 


3 


;i 


1)3 


2. 


(3,1,3) 


a 


= 14,/?. 


= -8 


A011922 


1, 3, 33, 451, 6273, 87363, 1216801, 16947843, . . . 


[1,3;4] 


4 


1 


ly 


3. 


(4,1,4) 


a 


= 23,/? = 


= —15 




1, 4, 76, 1729, 39676, 910804, 20908801, . . . 


[1,4;5] 


5 


'1 


1) 


4. 


(5, 1, 5) 


a 


OA R - 


- 0/1 




1, 5, 145, 4901, 166465, 5654885, 192099601, . . . 


[i,5,bj 


6 


^1 


1) 


5. 


(3,-1,3) 


a 


— — 


Q 

o 


A r\r\r\A (ifi 
AUUU4DD 


— 1, 3, 15, 35, 63, 99, 143, 195, 255, 323, 399, . . . 




2 


1 


3) 


6. 


(4,-1,4) 


a 


— I ,P — 


1 K 
it) 


AUoiU / o 


— 1, 4, 44, 319, 2204, 15124, 103679, 710644, . . . 




3 


'1 
^1 


/1^3 
4) 


7. 


(5,-1,5) 


a 


1/1 Q 
= 14, p = 


= z4 




— 1, 5, 95, 1349, 18815, 262085, 3650399 


[-1,5;4J 


4 


'1 
^1 


5j 


8. 


(6,-1,6) 


a 


= 2d,p = 


= OO 




— 1, 6, 174, 4031, 92574, 2125206, 48787199, . . . 


[-1,6;5J 


5 




6)5 


Bicentral sequences 


















9. 


(2,1,1,2) 


a 


= 4,/3 = 


-1 


A101265 


1, 2, 6, 21, 77, 286, 1066, 3977, 14841, . . . 


[i,i;3] 


2 


;i 


1)3 


10. 


(3,1,1,3) 


a 


= 6,/3 = 


-2 


A011900 


1, 3, 15, 85, 493, 2871, 16731, 97513, 568345, . . . 


[ia;4] 


2 


;i 


1)4 


11. 


(4,1,1,4) 


a 


= 8,/3 = 


-3 




1,4, 28, 217, 1705, 13420, 105652, 831793, . . . 


[i,i;5] 


2 


1 


1)5 


12. 


(5,1,1,5) 


a 


= 10,/? = 


= -4 


A054318 


1, 5, 45, 441, 4361, 43165, 427285, 4229681, . . . 


[i,i;6] 


2 


;i 


1)6 


13. 


(3,2,2,3) 


a 


= 3,/3 = 


-1 


A032908 


2, 3, 6, 14, 35, 90, 234, 611, 1598, 4182, 10947, . . . 


[2,3;1] 


5 


1 


2)1 
2)1 


14. 


(5,2,2,5) 


a- 


=5, /3= - 


3 




2, 5, 20, 92, 437, 2090, 10010, 47957, 229772, . . . 


[2,5;1] 


7 




15. 


(0,0,1,3) 


a 


= 2,/3 = 


1 


A000217 


0, 1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 66, . . . 


[0,1 ;4] 


1 


1 


1)4 


16. 


(0,0,1,4) 


a 


= 3,/3 = 


1 


A027941 


0, 1, 4, 12, 33, 88, 232, 609, 1596, 4180, . . . 


[0,1;5] 


1 


'1 


If 


17. 


(0,0,1,5) 


a 


= 4,/3 = 


1 


A061278 


0, 1, 5, 20, 76, 285, 1065, 3976, 14840, 55385, . . . 


[0,1;6] 


1 


;i 


If 



Table 4.1: Examples of symmetric lens sequences 





seed 


constants 




sequence 


label and symbol 




la. 


(3,1,2) 


a 


= 10, /3 = 


-- -5 


24, 3, 1, 2, 14, 133, 1311, 12972, 128404, 1271063,... 


[1,2;4] 


3(1,2)4 


b. 


(2,1,3) 








14, 2, 1, 3, 24, 232, 2291, 22673, 224434, 2221662,... 


[1,3;3] 


^(1,3)3 


2. 


(5,3,6) 


a 


= 6, /3 = 


-7 


. . . , 108, 20, 5, 3, 6, 26, 143, 825, 4800, 27968, ... 


[5,3;1] 


«(l,3)i 


3. 


(3,-1,4) 


a 


= 4, /3 = 


11 


. . . , 403, 104, 24, 3, -1, 4, 28, 119, 459, 1728,. . . 


[-1,4;2] 


3(1,4)2 


4. 


(15,12,20) 


a 


= 4, /3 = 


-13 


. . . , 400, 112, 35, 15, 12, 20, 55, 187, 680, 2520, ... 




2(4,5)3 


5. 


(21,6,10) 


a 


= 10, /3 = 


= -29 


. . . , 16796, 1700, 175, 21, 6, 10, 65, 611, 6016, . . . 




^(2,5)3 


6. 


(1,2,4) 


a 


= 5/2, (3 


= 


. . . , 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, . . . 


[l,2;3/2] 


3(1,2)2/ 



Table 4.2: Examples of non-symmetric sequences. Example 6 ( A000079 ) is integer in one 
direction only. 



seed constants sequence label and symbol 



1. 


(2, -1,2) 


a 


= -1, P 


= 3 


2, 2, -1, 2, 2, -1, 2, 2, -1, 2, 2, -1, 2, 2, . . . 


[-l,2;l] 


1(1,2)1 
1(1,2)2 
2(5,7)1 

1(1,1)2 

1(1,2)3 
3(5,2)1 


2. 


(3, -1, 2) 


a 


= 0, /3 = 


5 


2, 6, 3, -1, 2, 6, 3, -1, 2, 6, -1, 2, 6,-1, ... 


[-l,2;2] 


3. 


(14, -6, 15) 


a 


= 0, /3 = 


29 


14, -6, 15, 35, 14, -6, 15, 10, 35, 14, -6 . . . 




4. 


(1, 1, 0) 


a 


= 0, /3 = 


1 


1, 1, 0, 0, 1, 1, 0, 0, 1, 1,0, 0, 1, 1, 0, 0, 1, 1, ... 


[0,1;2] 


5. 


(4, -1, 2) 


a 


= 1, /3 = 


7 


2, 10, 15, 12, 4, -1, 2, 10, 15, 15, 12, 4, -1 . . . 


[-1,2;3] 


6. 


(10, -6, 33) 


a 


= 1, /? = 


49 


33, 88, 104, 65, 10, -6, 33, 88, 104, 65, 10, -6, ... 





Table 4.3: Examples of periodic sequences. Example 4 is known as A021913 . 
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Figure 4.2: Geometric representations (only general shape) for examples from Table 4.1: 
(a) Examples 1-4, (b) Examples 9-14, (c) Examples 5-8, (d) Examples 15-17. 




Figure 4.3: Geometric representations for examples from Table 4.2: (a) Examples 1-2 
(b) Example 3, (c) Example 4. 




Figure 4.4: Geometric representations for examples from Table 4.3: (a) Example 1, (b) 
Example 2, (c) Example 3, (d) Example 4. 





constants 




OEIS 


sequence 


label and symbol 




1. 


a — —3, 13 — 


1 


A001654 


0, 1, -2, 6, -15, 40, -104, 273, -714, 1870, -4895, 12816, . . . 


[0,1;-1] 


^(1,1)-^ 

-1(1,2)1 

2(i,i)-i 

^(l,3)-i 

-3(1,1)1 

-2(l,l)-2 

-4(1,1)1 


2. 


a — —3, (3 — 


5 


A075269 


2, -3, 12, -28, 77, -198, 522, -1363, 3572, -9348, . . . 


[2,2;1] 


3. 


a = -4, /? = 


3 




1, -2, 10, -35, 133, -494, 1846, -6887, 25705, -95930, . . . 




4. 


a = -4, (3^ 


-1 


A109437 


0, -1, 3, -12, 44, -165, 615, -2296, 8568, -31977, . . . 


[0-l;2] 


5. 


a = —5, (3 = 


1 


A099025 


0, 1, -4, 20, -95, 456, -2184, 10465, -50140, 240236, . . . 


[0,l;-3] 


6. 


a = —6, (3 = 


-4 


A084159 


1, -3, 21, -119, 697, -4059, 23661, -137903, 803761, . . . 


[i,i;-2] 


6. 


a = —6, (3 = 


1 


A084158 


0, 1, -5, 30, -174, 1015, -5915, 34476, -200940, . . . 


[0,1 ;-4] 



Table 4.4: Examples of formal lens sequences 
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But the above types of sequences do not exhaust the possibihties, as this example shows: 

. . . 2331, 407, 77, 21, 15, 35, 161, 897, 5187 . . . , 

which is a lens sequence with recurrence formula 6„ = Qb^-i — bn-2 — 34. We will arrive 
at a general rule that produces all integer lens sequences and an improved version of the 
integrality criterion in the last section. 



Invariants 

The formulas for the coefficients a and (3 in the recurrence formula (4.1) may be represented 
diagrammatically as shown in Figure 4.5, which exhibits the coefficients' algebraic "struc- 
ture" (the mnemonic role aside). The dots on the line represent the consecutive terms of the 
sequence. The arcs represent products of the joined terms, and the position above/below 
the line position indicates their appearance in the numerator/denominator of the formula. 
Dotted lines are to be taken with the negative sign. 



a+2 a p 




Figure 4.5: Diagrammatic representation of the lens sequence invariants 

Since the formulas do not depend on the particular choice of the three seed circles, one 
may position the diagram at any place in the line /sequence. In this sense, a and /3 represent 
invariants of the sequence with respect to translation along the sequence. But this also 
means that each of them gives rise to a new non-linear recurrence formula! 

Additional remarks on the Binet-like formula 

Denoting "jumps" around the central element bo by A+ = bi — bo and by A_ = bo — we 
get a more suggestive form of term w in the Binet-like formula for lens sequences, namely: 

A . - A_ A I + A_ 



w 



Va^ - 4 . 



2(a-2) a^-4 
Expressing a in terms of A we also obtain 

A, -A_ A, +A_ 



w 



2(A + A-i-2) VA2 + A~2 
Other representations of the formula for w include: 

bo , /3 A 6i - 6_i 



2 a-2j 2a^-8 

(a + 2)(6i -26o + 6_i) + (6i - 6_i)v/^^^ 

w = 

2(a2 - 4) 
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Note that re-indexing the sequence so that another term of becomes the central "60" 
will change the value of w in the formula (4.6). 

In order to better understand the situation, note the following simple general rule for 
this type of recurrences: 

Theorem 4.9. Let a sequence he given by the following Binet-like formula 

a;„ = acu" + + c (4.12) 

for some constants a, h, c and u). Then the sequence satisfies a non-homogeneous 3-term 
linear recurrence formula 

Xn+2 = {UJ + l/Lj)Xn+l - Xn + c{2 - U - l/u) . 

Moreover, if \uj\ > 1 then lim = uj, and for large n we have Xn — auj"^ + c. 
Proof. Direct. First note that (4.12) for the {n + 2)-nd term gives 

Adding these two gives, after some simple algebraic operations: 

Xn + Xn+2 = (c^ + l/uj)Xn+l + c{2 - UJ - I / Uj) 

which is equivalent to (4.12). □ 

Looking at the above, one hardly escapes the thought that lens sequences could be 
"explained" in terms of Chebyshev polynomials. This path did not however return any 
deeper insight. Instead, consider the following. 
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5 Underground sequences 



There are still some mysteries in the structure of lens sequences. One of the most remarkable 
properties is this: the entries of an integer lens sequence are products of consecutive pairs of 
a certain "underlying" integer sequence. Here is an example of a sequence of Vesica Piscis 
(A011922, see Example 1.1 in Section 1): 

b,: ... 3 1 3 33 451 6273 87363 1216801 ... 

AAA A A A A {["^^''ipiy 

Xi. 1 1 3 11 41 153 571 2131 



Thus the lens sequence may be represented as bi = fi-ifi, for some integer sequence (/j). 
Sequence (/j) will be called in this context the underground sequence of sequence 
(Table 6.1 provides examples). We present the formalism of this amazing and unexpected 
property. 

Let us start with a general fact. 

Theorem 5.1. Any factorization {/«} of lens sequence in a sense that bn = fn-ifn satisfies 
3-term recurrence formula 

fn+2 + fn-2 = Oifn (5.1) 

Proof. Starting with the expression for a given in Proposition 3.1, we have 

^ _ bn-l ^ bn+2 _ fn-2fn-l ^ fn+lfn+2 _ fn-2 ^ fn+2 _ fn-2 + fn+2 
bn fn—lfn fnfn+1 fn fn fn 

□ 

It should be borne in mind that this property is true for any — not necessarily integer — 
factorization of {bi}. Such factorizations are easy to produce, e.g., set /o = 1, /i = bi, 
f2 = &2/&1, /a = b^bi/bo, fi = etc. However: 

Theorem 5.2 (Factorization theorem). Any integer lens sequence {bn) may be factored into 
an integer sequence (fn) so that bn = fn-ifn- If the lens sequence is primitive, the factoriza- 
tion is — up to a sign — unique. Moreover, in such a case = gcd(6„, ^n+i)- 

Proof. Assume that {bn) is a primitive lens sequence. Consider three consecutive terms 
(a, 6, c) and define 

/ a f u\ f ^ ^ f ^ gcd(a,6)c 

= A( h\ = Scd a, b) f2 = — = /3 = — = 

gcd(a,6) /i gcd(a,&) /s b 

Clearly, a = /o/i, b = /1/2, and c = 72/3. We need to show that these four terms are 
integers. Terms /o, /i, and /2 are integer by definition. As to the last term, use the formula 

b'^—ac . 



^ o-c ^ gcd(a6)c ^ 

b b 
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hence /s is an integer. Thus /o, /i, /2, /s G Z and the integrahty of the whole sequence (/„) 
follows immediately from 5.1. 

As to uniqueness of factorization of a primitive lens sequence, assume a contrario that 
two integer quadruples, {/o, /i, /2, fs} and {go, gi, g2, gs} are the initial terms of two different 
factorizations of Then go/ fo = p/q for some mutually prime p,q & At least one of 
p and q is not 1; assume that it is p 7^ 1. Since gigi+i = fifi+i, we must have 

{90, 9i, 92, 93} = {- fo, - /i, - /2, - /s} C Z . 
q p q p 

Thus q\fo and gl/2 (because gcd(p, g) = 1). But this means that q \ a (since a = /o/i), q \ b 
(since 6 = /i/2), and q \ c (since c = /2/3), against the assumption of primitivity of the lens 
sequence. □ 

The 3-term recurrence 5.1 for the underground sequence involves only a. Another inter- 
esting non-linear 4-term recurrence involves only j3: 

Proposition 5.3. Any underground sequence {/j} of a lens sequence {bi} satisfies the fol- 
lowing quadratic recurrence formula: 



det 

Proof. For any n we have 



fn fn+1 
/n+2 fn+3 



fn+3fn — fn+lfn+2 — — (5.2) 



r p n n fn fn+1 f n+2 fn+3 n r 

JnJn+3 ~ Jn+lJn+2 — 7 7 Jn+lJn+2 

Jn+lJn+2 

bn+lbn+3 , bn+lbn+3 — ^^+2 n 

= -T On+2 = 7 = -P 

On+2 0„+2 



□ 



Note that not any initial quadruple (/i, /2, /s, /4) leads via recurrence 5.1 or 5.6to a 
sequence that underlines a lens sequence. When do they? First, we notice that the under- 
ground sequences of lens sequences have an interesting anatomy. It turns out that they are 
determined by three-term linear recurrences with variable "constants" . Here is the central 
theorem for the underground sequences: 

Theorem 5.4 (Underground Sequence Structure), (i) Let k,s E Ij be two constants. Define 
a sequence f by 

J ^ f fc/n-i - fn-2 if n is even 
\sfn-i - fn-2 if n is odd. 

with some arbitrary initial terms fo, fi G Z. Define bn = fn-ifn- Then (6„) is a lens 
sequence. The constants of its recurrence formula 

bn = abn-i -bn-2 + P (5.4) 
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are 



a = ks — 2 

f3 = kff + sfl~ksfoh 



(5.5) 



(a) Every lens sequence is of such type. In particular, for a primitive lens sequence with a 
seed (6-1, 605^1) = {(i,b,c) the underground sequence is defined 



fo = gcd(a,6), 
a + b 



fi 



2 ' 



/i = gcd(6,c), 
^" fl • 



Proof. We start with part (ii). Let (/j) be a sequence defined by 5.3. First, we shall show 
that the following expression 



A, 



det 



fn—3 fn—2 
fn— 1 fn 



(5.6) 



is an invariant of a sequence (5.3), that is it does not depend on n. Indeed, let p denote k 
or s, depending on whether n is even or odd (it will not matter!). Then 



An = det 


fn—3 
fn-1 


fn-2 
fn 


= det 


fn—3 
fn-1 


Pfn-3 
Pfn-1 


fn 
fn 


= det 


fn—3 
fn-1 


fn-4 
~fn-2 


= det 


fn-2 
fn-4: 


fn-1 
fn—3 





A 



n-l 



Thus, by induction, the value of A„ does not depend on n. Hence it may be brought down 
to the first two terms of the sequence and, after simple substitution, shown to be 

An = kf^ + sf^-ksfofl. 

As to the recurrence formula, calculate the sum + bn+i '■ 

bn+l + bn-l = fnfn+l + fn-2 fn-1 = {o-fn-1 ~ fn-2)(bfn — fn-l) + fn-2 fn-1 
= ks fn-lfn — kfn-i — sfnfn-2 + fn-lfn-2 + fn-lfn-2 
= ks fn-lfn — fn-l{kfn-l ~ fn-2) — fn-2{sfn ~ fn-l) 

= ks fn-lfn — fn-lfn — fn-2fn+l (5.7) 
= ks fn-lfn — 2fn-lfn + fn-lfn ~ fn-2fn+l 
= {ks — 2) fn-lfn + {fn-lfn — fn-2fn+l) 
= {ks - 2) bn- An. 

Now, rename ks — 2 = a, and A„ = — /3 (as it does not depend on n). Then (5.7) becomes 
an inhomogeneous three-term recurrence formula 

bn+i + bn-1 = abn + p. 
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Finally, we need to show that is actually a lens system sequence. As for a, we need 
simply to show 



a 



bn-i _^ bn+2 

bn bn+1 



(see Proposition 3.1). Consider the right hand side: 



bn-i _^ bn+2 



fn-2fn-l /n+l/n+2 _ fn-2 _|_ fn+2 fn-2 + fn+2 



bn bn+l fn—lfn fnfn+1 fn fn fn 

where the last step is true because of (5.3). Indeed, assume n is even: 

fn-2 + fn+2 fn-2 + ^/n+l — fn 



a , 



fn 



fn 

fn-2 + k{sfn — fn-l) ~ fn 
fn 

ksfn kfn-l fn-2 fn 
fn 

ksfn fn fn 



fn 



ks 



a . 



As to the formula for j3 in terms of a seed, follow similar calculations as in the proof of 
Proposition 5.3: 

^ bib, - bl 

12 = —. = -, = -P, 



A, = det 



/o /i 
/2 /a 



/1/2 



This ends the proof of (i). The proof of (ii) follows easily. Let (a, 6, c) be a primitive seed of an 
integer sequence. Define a quadruple of numbers (/o, /i, /2, /a, fi) by setting /i = gcd(a, 6), 
/2 = gcd(5,c), and /o = a//i, and f^ = C//2. Then k = (/o + /2)//i and s = Ui + f^)/f2- □ 



Remark 5.5. Note that Aj may be viewed as a quadratic form given by the matrix 

G- 



k ks 
s 



evaluated on the vector v = [/o,/i]"^, i.e., A„ = v-^Gv. In particular, vectors (/j, /j+i) G M 
all stay on a quadratic defined by G. 

As a corollary of the above theorem, we arrive at a simple criterion on whether a triple 
of integers is a good candidate for an integer lens sequence, improving that of Proposition 
4.2: 

Theorem 5.6 (Integrality Criterion 2). A triple {a,b,c) G Tj is a seed of an integer lens 
sequence iff 

(z) gcd(a, b) gcd(6, c) = 6gcd(a, b, c) 
(ii) {gcd{a,b)Y divides (a + 6) gcd(a, 6, c) 
(Hi) (gcd(6, c))^ divides (6 + c) gcd(a, 6, c) 
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a+2 

ir\nr\\ / n \ 




Figure 5.1: Diagrammatic representation of the lens sequence invariants, calculated from the 
underground sequence (cf. Fig. 3.2) 

The above result leads to another property of lens sequences: 

Corollary 5.7. The sum of any pair of consecutive terms of a lens sequence is a multiple 
of a square, namely: 



bn + bn+i 



^fn ^/ ^ even _ \ k ■ square if n is even 
sf^ if n is odd \ s ■ square if n is odd. 



Proof Elementary: bn + 6„+i = fnfn-i + fn+ifn = fnifn+i + fn-i) = fnPfn, where p stands 
for k or s, depending on the parity of n. □ 

For example, the sums of two consecutive entries of A011922 are perfect squares: 

... 3 1 3 33 451 6273 87363 1216801... 

YVVYV Y Y R.* 

4 4 36 484 6724 93636 1304164 v 
2^ 2^ 6^ 22^ 82- 306^ 1142^ 

The sequence A101265 (of label [1,1; 3]) gives: 

... 1 1 2 6 21 77 286 1066 3977 



YYYYYY Y Y {1 

2 3 8 27 98 363 1352 5043 
2-1^ 3-1- 2-2^ 3-3^ 2-7^ 3-11^ 2-26^ 3-41- 



add 



Generating lens sequences 

Let us return to the question of generating integer lens sequences. In Section 2, we considered 
triplets of the form (a, b; k) that label a large family of lens sequences (see (3.1)), but not all 
of them. The existence of underground sequences allows one to label all lens sequences. 

Definition 5.8. A symbol of a lens sequence is the quadruple ^{p,q)^ which defines the 
underground sequence (/«) with /o = p, fi = q, and with constant s and k in (5.3), and 
therefore defines the corresponding lens sequence (6^), namely, 6j = fi^ifi- More directly, 
symbol ''(p, qY defines a lens sequence via its seed (a, b, c) = {{sp — q)p,pq, q{kq — p)). 
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Any integer quadruple ^{p, leads to an integer lens sequence. And vice versa, given a 
seed of a primitive sequence (a,6, c), we easily reproduce the symbol: 



p = gcd(a, b), k 
g = gcd(6,c), s 



b + c 

a + b 
p2 



Proposition 5.9. The lens sequence generated by q) is primitive if and only if 

gcd(p, q) = gcd(p, k) = gcd(s, q) = 1. 



(5i 



Proof. Write the "central" four terms of the underground sequences and the corresponding 
lens sequence: 

ifi) ■ •••,/-! = {sp ~ q)Jo = P, fi = qj2 = {kq - p), . . . 
(bi) : ...,a = {sp- q)p, b = pq,c = q{kq - p), . . . 

For to be primitive we must have gcd(a, 6, c) = 1, which implies (5.8). □ 

Corollary 5.10. If the underground sequence (fi) contains p = ±1, then the corresponding 
lens sequence (bi) admits label {a,b;k). More precisely: 

(a, b; k) corresponds to the symbol """"^(l, b)^ ; 
the symbol ^{1, b)^ corresponds to (s — b, b; k) . 

Remark 5.11. To use this generator of sequences as a unique label system for lens sequences, 
one would have to remove the ambiguity of the choice of the initial terms. We may demand 
that, say, p = fo has the smallest absolute value among (/j) and that > /o ^ fi- 

Remark on diagrammatic use of symbols. The first of the following two diagrams 
means that 13 is obtained as 13 = 5 x 3 — 2. The second represents equation 2 = 5x3 — 2: 



13 



13 



Now, the symbol ^(2,3)^ may be graphically developed into an underground sequence (/j 
and consequently into a lens sequence in the following way: 



start with these two 



^1 ^ 




multiply consecutive pairs 
42 7 -2 6 39 130 370 



(Xi) 



(bd 
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Note that the four central terms of (/,) suffice to generate the three central terms of 
which yield the constants a and (5. The recurrence formula for sequence (/j) is bilateral and 
may be represented diagrammatically as shown below: 




... -6 



-7 



-12 3 1 



13 



10 



37 ... 




Example 1 (Vesica Piscis): b = A011922 f = A001835 = A079935 
{6,} -(..., 1, 3, 33, 451, 6273, 87363, . . .), &„+i = 146„ - 6„_i - 8 

{/i} = (••■, 1, 3, 11, 41, 153, 571, . . . ), ^(1, 1)^ (number of domino packings in a (3 x 2n) rectangle) 

Example 2 (Golden Vesica): b = A064170 f = A001519 

{&J = (..., 1, 2, 10, 65, 442, 30, 26, 20737, . . .) &„+i = 76„ - 6„_i - 3 

{/J = (...,1,2,5,13,34,89,233,...), ^(1,1)^ (odd Fibonacci numbers) 

Example 3: b = A000466 f = A005408 

{h} = (..., -1,3, 15, 35, 63, 99, 143, 195, 255, . . .) 6„+i = 2b„ - fe„_i + 8 
{/J = (.••- 1,1,3,5,7,9,11,13,15,...), 2(-l,l)2 (odd numbers) 

Example 4: b = A081078 f = A002878 

{6,}- (...,-1,4,44,319,2204,...) = 76„ - 6„_i + 15 

{/,} ^ (... - 1,1,4,11,29,76,199,521,1364...) , 3(_i^i)3 (^^j^ L^^as numbers) 
Example 5: b = A005247 f = A005247 

{&J = (...,2, 2, 3, 6, 14, 35, 90, 234, 611, 1598, . . .) 6„+i = 3fe„ - 6„_i - 1 

{/J = (..., 2, 1, 3, 2, 7, 5, 18, 13, 47, 34, ...),, ^ ^(2, 1)^ 

(even Lucas numbers interlaced with odd Fibonacci numbers) 

Example 6: b = A027941 f = A005013 

{&,} = (...,0, 1, 4, 12, 33, 88, 232, 609, . . .) 6„+i = 36„ - 6„_i + 1 



Table 5.1: Examples of underground sequences f for some integer lens sequences b. 

Towards the meaning of the underground sequence 

Finally, let us consider yet another recurrence formula for lens sequences and its surprising 
context. Let us start with this: 

Proposition 5.12. Lens sequences obey the following nonlinear 4-step recurrence formula: 




Example 7: b =A000217 f = A026741 

{&,} = (...,0, 1, 3, 6, 10, 15, 21, 28 . . .) 6„+i = 26„ - b^-i + 1 

{/,} = (...,0,1, 1,3,2,5,3,7,4,9,...) , \0,1)^ 



(the sequence of natural numbers interlaced with odd natural numbers) : /., 




'n+2 — 



{hn+l-P){hn-f3) 
bn-1 



(5.9) 
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Proof. Rewrite the definition of (3 in tlie form ac = b{b — (3). Since the three consecutive 
terms (a, b, c) may start with any entry of the sequence, let us write its two instances: 

bn+lbn-l = bn{bn - 
bn+2bn = bn+l{bn+l — P) 

Muhiply side-wise to get &n+i&n-i&n.+2&n = — /5)(fen+i — (3)bnbn+i- Cancehng the repeated 
terms results in 

bn-lbn+2 = {bn " (3){bn+l " 

which is equivalent to (5.9). □ 
Corollary 5.13. A lens sequence satisfies the following identity: 



det 



bn+l ~ 13 bn+2 
bn~l bn - 13 



0. 



This leads to yet another implication. The above determinant (5.9) may be written in 
the form: 

bn+l bn+2 
bn-1 bn 



det 







" 














which looks like a characteristic equation with (3 playing the role of the eigenvalue. Note 
that it does not depend on n. What are the corresponding eigenvectors? 

Theorem 5.14. If (fi) is the underground sequence of a lens sequence (bi), then the following 
"eig en- equation" holds: 



bn+l bn+2 
bn-1 bn 



fn+1 
^fn-1 



/3 



fn+1 
~fn-l 



Proof. Write 6's in terms of /'s and use the definition of (3 written in terms of /'s. 



□ 



6 Summary 

A lens sequence is an integer sequence (bi) that satisfies two conditions: 

(i) bn = Oibn~i — bn~2 + [3 [recurreucc formula] 

(ii) + b^ = a ab + [3{a + b) [compatibility relation] 

where a and (3 are constants and a and b are any two consecutive terms of (bi). These two 
conditions assure that the sequence has a geometric realization in terms of the curvatures of 
a chain of circles inscribed in a symmetric lens (the space of the overlap of the interiors or 
exteriors of two congruent circles). The sequence constants may be viewed as invariants of 
a process i ^ b^. They may be calculated from a seed, i.e., any three consecutive sequence 
terms (a, 6, c): 

ab + bc + ca , ^ — ac 
a = 1 and p = — ; 
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or, for a, alternatively 



a — — h 



Other, nonlinear, recurrence formulas for the lens sequence include: 

[two-step formula] 26„+i = 6„q; + /3 ± ^J{a^ - 4) 6^ + 2{a + 2)/3 6„ + (5"^ 

[three-step formula] hn+ihn + bn+ibn-i + bnbn-i ^ {a + 1)6^ (only a) 

[three-step formula] bnbn — &n+i^n-i = Pbn (only (3) 

[four-step formula] bn+2bn-i = (&n+i - /3)(&n - 

[four-step formula] -|- bnbn-2 = ab^bn-i 

The sequence constants have a geometric meaning: a codes the angle under which the 
circles forming the lens intersect (if they do), or, more generally, the Pedoe product of the 
lens circles. The value of 13 reflects the size of the system. There are two basic properties 
determined by geometry: (a) the sum of the inverses is determined by the length of the lens, 
and (b) the limit of the ratio of consecutive terms is determined by the aforementioned lens 
angle: 



E2 Voi^ - 4 , bi+i a + \fa'--^ 
— — — and lim — — = = A , 
bn -(5 i^oo hi 2 



The number A, the characteristic constant of {hi), allows one to express the lens sequence 
by a Binet-type formula 

bn = wA" + wX" + 7 

where 

a — 2h-\- c c — a i-^ — /3 

w^— — ^ — VQ!"'-4, 7 



2(q; - 2) 2(q;2 - 4) ' ' a - 2 ' 



and where the bar denotes natural conjugation in the field Q(-\/ci;^ — 4) . The constant A 
is an example of a (quadratic) Pisot number, an algebraic integer, the powers of which 
approximate natural numbers. In particular, 6„ fa u'A"'-|-7. Lens sequences can be expressed 
also as combinations of Chcbyshev polynomials. 

The most mysterious property of a lens sequence is that its terms may be formed by taking 
products of pairs of consecutive terms of another sequence. This "underground" sequence 
has an alternating recurrence rule, different for odd and even terms. Namely 6„ = fn-ifn , 
where: 

kfn-i - fn^2 if n is even 
sfn^i - fn-2 if n is odd. 

The constants of the sequence (6^) may now be expressed as 

a — ks — 2 

P-sf^ + kfi-ksfoh. 



fn 



28 



It follows that the integer lens sequences may be determined by four arbitrary integers 



^(/o,/i)^ 



Choosing for /q the term with smallest absolute value allows one to treat the above quadruple 
as a symbol that labels the corresponding lens sequence. 

The underground sequences automatically satisfy the following two recurrence formulas: 



(i) fn+2 + fn-2 = a/n 
{ii) fn+2fn-l — fnfn+1 



-/3 



or 



det 



fn— 1 fn 
fn+1 fn+2 



More precisely, the set of the sequences that are underground sequences for lens sequences 
coincides with the intersection A n A, where A denotes the space of sequences satisfying 
linear recurrence (i), and A denotes the set of sequences satisfying (ii). 

An intriguing property holds — the eigenvectors of matrices assembled from the terms 
of a lens sequence are vectors with entries from the corresponding underground sequence: 



bn+l bn+2 
bn-1 bn 



fn+1 
— fn-1 



fn+1 
— fn-l 



But the full meaning of the underground sequences remains to be understood. 



Appendix 

Below, we summarize the general formulas for each of the five types of symmetric integer 
lens sequences. Recall that L — length of the lens, R — radius of the lens circles, 5 — their 
relative distance, and A = characteristic constant. Due to symmetry, the sums of the recip- 
rocals are curtailed to one (right) tail of the sequence. 



1. Seed: [n, l,n] . Symbol: ""+^1, 1)"+^ Recurrence: a = (n+iy-2, (3 = l-n^ 
Geometry: L = 2^, (5=^. (Inner chain) 

Characteristic constant: A = (^iWpi^y ^ n^ +2n-i+in+i)./i^^^^ 

Binet: bk = ^"+^^+"+1 

Sum of reciprocals: ^ l/bk = ^ + i + ■ ■ ■ = \ + \J ^ 
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2. Seed: [n, 1, l,n] . Symbol: ^(1, 1)"'+^. Recurrence: a — 2n, (3 — 1 — n . 



Geometry: i2 = 2^, L = 2^^, 5 = -^=4=. (Inner chain) 
Characteristic constant: A = (^ Vw+^V^ y = n+v|!Ei 

OO / 

Sum of reciprocals: ^ 1/6^. = 1 + i + . . . = . 

! 1 V 



3. Seed: [n, 2, 2, n] . Symbol: ^(2, 1)"+^. Recurrence: a^n , P^2-n . 
Geometry: = L = 2v^ = 2^, 5 = y=. (Inner chain) 



Characteristic constant: A = 
Binet: = i (l + + i (l - ^if) A-'^ + 1 

OO I 

Sum of reciprocals: E V^fe = i + ^ + ■ ■ ■ = Iv ^ 

fc=i V 



4. Seed: [n, -l,n] . Symbol: "-^(1, 1)"+^ Recurrence: a = (n-lf-2, (3 = n^-l. 
Geometry: L = 2^, 5=;^. (Outer chain) 

Characteristic constant: A = (^!!zi±Vp(^^ ' ^ 

Binet: 6. ^ '^^IT"' 

OO 

Sum of reciprocals: E = n + ■ ■ ■ ^ 



5. Seed: [0, l,n] . Symbol: ^(1, 1)"+^ Recurrence: a^n-l , = 

Geometry: R——{n+l), L = 2a/ (n + l)(n — 3), 5 = ;^7^3=- (Outer chain) 

/ n+l+^(n-3) V _ n-l+^(n+l)(n-3) 



Characteristic constant: A 



Binet: 6^ = i ( ^ + , ^ A*^ + | ( ^ - , " 1 A-'^ + -\ 

" 2 yn-3 y'(n+l)(n-3) y 2 \^n-3 ^(„+l)(„-3) ^ "-3 

Sum Of reciprocals: £ 1/6^ = 1 + i + . . . = ^+i-V(^+^)("-3) 

fc=0 " 
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